Sensitive imaging of electromagnetic fields with paramagnetic polar molecules 



Sergey V. Alyabyshev 
Department of Chemistry, University of British Columbia, Vancouver, BC V6T IZl, Canad^ 

Mikhail Lemeshko 
ITAMP, Harvard-Smithsonian Center for Astrophysics and Harvard 
Physics Department, Cambridge, Massachusetts 02138, USA 

Roman V. Krems 

Department of Chemistry, University of British Columbia, Vancouver, BC V6T IZl, Canada^ 

(Dated: February 10, 2012) 

We propose a method for sensitive parallel detection of low-frequency electromagnetic fields based 
on the fine structure interactions in paramagnetic polar molecules. Compared to the recently imple- 
mented scheme employing ultracold *^Rb atoms [Bohi et al., Appl. Phys. Lett. 97, 051101 (2010)], 
the technique based on molecules offers a 100-fold higher sensitivity, the possibility to measure both 
the electric and magnetic field components, and a probe of a wide range of frequencies from the dc 
limit to the THz regime. 



Sensitive detection of weak electromagnetic fields is 
critical for many applications ranging from fundamen- 
tal physics measurements [T] to biomagnetic imaging of 
the brain and heart to the detection of explosive ma- 
terials [3j. A tremendous progress in measuring mag- 
netic fields has been recently achieved leading to the de- 
velopment of Hall effect sensors [4 , SQUID sensors [5] , 
force sensors [Sj, sensors based on microelectromechani- 
cal systems "7; , and NV centers in diamond [5| , as well as 
atomic magnetometers p], making it possible to achieve 
the magnetic field sensitivity of 0.1 fT Hz^^^^ and to 
detect the magnetic field of a single electron, with steps 
being taken towards the detection of the magnetic field of 
a single nuclear spin [TU]. At the same time, the develop- 
ment of scanning capacitance microscopy [XT', scanning 
Kelvin probe [12] , and electric field-sensitive atomic force 
microscopy [13] advanced the techniques for measuring 
electric fields to the level of probing individual charges. 
An unprecedented accuracy of 10~^ electron charge was 
achieved with the use of single-electron transistors [H]. 
Yet, even with the sensitivity pushed to its fundamental 
limit, none of these methods allows for parallel measure- 
ments, i.e. imaging of the field amplitudes and phases at 
many spatial points at the same time. Recently Bohi et 
al. proposed to use ultracold atoms for sensitive paral- 
lel imaging of microwave fields with frequencies in the 
range 2.5 — 14 GHz [T^. The method relies on measur- 
ing the phase difference between two hyperfine states of 
®''Rb, accumulated due to an interaction with the mag- 
netic component of the microwave field. The acquired 
phase difference is proportional to the evolution time, 
the magnetic moment of the atom, and the magnetic field 
amplitude. However, long measurement times lead to im- 
age blurring due to the atomic motion and decoherence, 
resulting in a compromise between field sensitivity and 
spatial resolution. 

Although measuring the electric component of an ac 



field is several orders of magnitude more efficient than 
detecting the magnetic component [16] , atoms possess no 
permanent electric dipole moments, rendering the detec- 
tion of the Zeeman shift the only possible option. Here we 
describe a technique for parallel, noninvasive, and com- 
plete (amplitudes and phases) imaging of electromagnetic 
fields with an ensemble of polar open-shell molecules, 
many of which have been successfully cooled and trapped 
in experiments [TTl - EU] . We show that the presence of 
permanent electric and magnetic dipole moments and the 
variety of molecular rotational constants allow for the de- 
tection of both electric and magnetic field components, in 
a wide range of frequencies from the dc limit through the 
radio and microwave to the THz frequency range. We 
show that measuring the electric component of an os- 
cillating field must result in shorter measurement times 
compared to the atomic experiments, and consequently 
higher spatial resolution, which is mainly limited by the 
optical detection scheme and photon scattering. 

First, consider molecules with a dipole moment 
d, subject to known dc magnetic and electric fields, Bq 
and Eg, both pointing along the laboratory Z axis. The 
molecules are described by the following Hamiltonian: 

F = 6,N2+7N-S+5sMBBo-S-d-Eo, (1) 

where N and S are the rotational and spin angular mo- 
menta of the molecule, br and 7 are the rotational and 
spin-rotation interaction constants, /is is the Bohr mag- 
neton, and gs = 2.0023. In the absence of external fields, 
the states of a molecule are labeled by \N,J,M), 
where J = N -I- S is the total angular momentum and 
M is the projection of J on the Z axis. We use the 
same quantum numbers to label the molecular states in 
the presence of dc fields, bearing in mind that J is not 
conserved. 

The effect of combined Bq and Eq fields on the ro- 
tational states of a molecule is shown in Fig. [I] (a) 
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FIG. 1. (a), (b) Energy levels of the SrF(^E+) molecule 
{K = 7.53 GHz, 7 = 74.7 MHz) in an electric field of Eq = 
10 kV/cm as a function of magnetic field Bo; (c) Frequency 
dependence of the ac field sensitivity for SrF in a linearly 
polarized mw field for different electric fields; the red and 
black lines correspond to the 2—^3 and 1—^4 transitions 
respectively. The dashed line represents the sensitivity to the 
magnetic field component of the ac field that can be achieved 
in experiments with atoms [15j : (d) Same as in (c) but for the 
CaH(2E+) molecule (6^ = 128.3 GHz, 7 = 1.24 GHz) 



and (b) for the case of SrF (2^+) pD]. We assume 
that the molecules are initially prepared in the magnetic 
low-field seeking Af -component of the rotational ground 
state, state |2) = |0, 1/2,1/2) in Fig. [ija). This can be 
achieved by cooling molecules to subKclvin temperatures 
and confining a molecular cloud in a magnetic trap |17j . 
Alternatively, molecules can be confined in an electric or 
optical trap and transferred to state |2) by a sequence 
of microwave pulses [5T]. Molecules can be cooled by 
a variety of recently developed experimental techniques 
such as buffer-gas cooling. Stark deceleration, or laser 
cooling [HHID]. 

State 1 2) exhibits an avoided crossing with the mag- 
netic high-field seeking state |3) = |1, 1/2, 1/2) at the 
magnetic field B* = 659.5 mT. The states |2) and |3) 
have the opposite parity and, due to the spin-rotation 
interaction, represent linear combinations of the states 
with spin projections Mg = ±1/2 22 . Therefore the 
1 2) — > 1 3) transition is dipole-allowed and can be used 
to detect the electric component of a resonant if or mi- 
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FIG. 2. Sensitivity of the |2) -i> |3) transition in SrF(2S+) 
to electromagnetic fields as a function of the dc electric field 
Eo and the frequency of the microwave field. The peaks are 
labeled by the frequency of the ac field in MHz. The dot- 
ted line shows the sensitivity of the |2) — >■ |3) transition in 
SrF(^E"'") to the magnetic field component of the microwave 
field. 

crowave field. 

The measurement requires placing the molecular en- 
semble close to the source of the field to be measured, 
E{r,t) = l/2[E(r)exp(-iwi) + E*(r)exp(iwt)], and an 
addition of background dc magnetic and electric fields, 
with magnitudes Bq and Eg. By adjusting Bq and Eq, 
the energy splitting between the states |2) and |3) can be 
tuned in resonance with the field frequency uj. The trap- 
ping field, if any, must be switched off in order to allow 
the field E{r, t) to drive resonant oscillations between the 
states |2) and |3) during the free evolution time r. The 
Rabi frequency of the oscillations is given by 

n^{Bo, Eo) = E^{r)d^{Bo, Eo)/h, (2) 

where a = {+; 0; — } denotes the polarization of the oscil- 
lating field with respect to the laboratory Z axis, Ea{r) 
is the corresponding component of the field E{r,t), and 
da{BQ, Eq) is the transition dipole moment between the 
two states. In order to compute da{Bo, Eq), we expand 
the molecular states |n) plotted in Fig. [l] as follows: 

E 4%Ms(Bo,Eo)\NM^}\SMs}, (3) 

N,Mn;Ms 

which gives 

W (R \ - \^ \^ (2)* (3) 

N,Mn,Ms N',M'j^,M'g 

{NMn\Y,^\N'M'j^)Sms,m's- (4) 

Here, Mn and Ms denote the projections of the angular 
momenta N and S on the Z axis, respectively. 
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FIG. 3. The ac electric field sensitivity for SrF( E ) as a 
function of the spatial resolution, for n = 10^" cm~^ (solid 
lines) and n — 10^^ cm"'' (dashed lines). Different colors 
correspond to different ac field frequencies. 



After time r the probability to detect a molecule in 
the state |3) at the spatial point r is 



"3(r) 



n2(r) + (r) 



f]„(r)T 



(5) 



where n2 and are the densities of the molecules in the 
states |2) and |3). From the measured value of P3(r) one 
can calculate ^lai'r) and, using Eq. ([2|, the components 
of the electric field. The amplitudes, Ex{r), Ey{r), and 
Ez{r), and the relative phases can be reconstructed by 
measuring (r) with the background magnetic field Bq 
pointing along the a;, y, and z axes, by analogy with ref. 
|15j . From E(r), the spatial distribution of the magnetic 
field can be calculated using the Maxwell equations, while 
the amplitudes of the two fields are related as i? = cB. 

The single-shot sensitivity of the measurement to the 
ac electric fields is given by [IS]: 



Ve 



[V/cm ^ 



100 dM,Eo)V^y 



(6) 



where n is tlie density of molecules, Vcs — 27rcr^gp is 
the effective imaging volume, (Toff is the dispersion of the 
spatial coordinate, and p represents the 1/e radius of the 
cloud. 

The dependence of the transition dipole moment da 
on the background fields Bq and Eq renders the sensi- 
tivity frequency dependent, rj^ = •q'^iuj), as shown in 
Fig. fl] (c) by red curves for the |2) — > |3) transition in 
SrF(^+). One can see that the frequency dependence 
exhibits sharp minima featuring the sensitivities on the 
order of 10^^ — 10^'' V/cm Hz~^/^. The positions of the 
minima can be controlled by tuning the splitting of the 
levels 1 2) and |3) with the background electrostatic field 



Eq, and thereby shifted towards smaller frequencies |23j . 
The accessible frequency range can be extended by ini- 
tially preparing the SrF molecules in the high-field seek- 
ing state, |1) = |0, 1/2,-1/2), and driving the transition 
to the state |4) = |1, 3/2, —3/2). The frequency depen- 
dence of the sensitivity corresponding to the |1) -> |4) 
transition is shown in Fig. [l] (c) by black lines. Choosing 
a molecule with a different rotational constant allows for 
the detection of a completely different range of accessi- 
ble frequencies. As an example, the rotational constant 
of CaH(^E+) molecule [T^ is about 17 times larger than 
that of SrF, which gives access to microwave fields of 
frequencies w ~ 100 — 500 GHz, as shown in Fig. [T] (d). 

Interestingly, the transition dipole moment da{BQ, Eq) 
for the 1 2) — >■ |3) transition in molecules vanishes at 
certain combinations of Eq and Bq. At these particular 
combinations, the molecules become transparent to the 
resonant microwave field. This is demonstrated in Fig- 
ure [2] The figure illustrates that the magnitude of the 
dc electric field Eq, for which da{Bo,EQ) ~ 0, depends 
sensitively on the frequency of the resonant transition 
(which can be tuned by varying Bq). This can be used 
for sensitive detection of the magnitude Eq of a dc elec- 
tric field, given the magnitude of Bq and the frequency of 
the microwave field. Conversely, this can be also used for 
sensitive detection of the magnitude Bq of the dc mag- 
netic field, given the magnitude of Eq and the resonant 
microwave frequency. 

Longer interaction time t results in increased sen- 
sitivity to electric fields. The sensitivity is, however, 
gained at the expense of the spatial resolution that de- 
creases with r due to the molecular motion and de- 
coherence. The effective spatial resolution, SeS = 
1/2 

2(cr^ -I- cr?jjjg -f (Tpg) , can be calculated from the dis- 
placements (Tr =^ Tyj2kj^TJrn during the measurement 
time (t), (Tiing — T\^^^ ^ksT I m during the imaging pulse 
(Timg), and dps — frec'''img\/2rT\^/3 duc to the photon 
scattering. Here /cs is the Boltzmann constant, T is the 
temperature of the gas, m is the mass of a molecule, Urcc 
is the recoil velocity, and F is the scattering rate (for SrF, 
F = 3 MHz 20 ) . Fig. [s] illustrates the relation between 
the sensitivity and the spatial resolution for an ensem- 
ble of SrF molecules and different frequencies of the field 
detected. 

The energy level structure of paramagnetic molecules 
can also be used to probe sensitively static or off-resonant 
rf and microwave fields. This can be achieved by measur- 
ing the phase accummulation in a Ramsey-type sequence 
consisting of two 7r/2 pulses The first 7r/2 pulse pre- 
pares the molecules in the equal superposition of states 
1 2) and |3), which acquire a relative phase proportional 
to the Stark shift, A0 — dcs{BQ, Eq)Et, due to the ef- 
fective dipole moment doff (^Oi -E'o) = '^^\iE^'' during the 
evolution time r. The second 7r/2-pulse transforms the 
relative phase into a population difference. Because the 
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FIG. 4. (a), (b) Energy levels of the NH(3E-) molecule [h^ = 
490.0 GHz, 7 = -1.65 GHz, and A = 27.6 GHz) as a function 
of the magnetic field Bo\ (c) Frequency dependence of the ac 
field sensitivity. The background electric field Eq — 0. 



d ^ 1 a.u. amounts to 
the proposed scheme. 

In summary, we have described a technique for sensi- 
tive parallel measurements of electric and magnetic field 
components of electromagnetic fields, both dc and oscil- 
lating with frequencies ranging from a fraction of a kHz to 
THz. The method, based on tunable energy level struc- 
ture of paramagnetic molecules in superimposed electric 
and magnetic fields, allows one to achieve the sensitivity 
on the order of ^V/cm Hz^i/^ ^nd 100 fT Hz^i/^ foj. ^he 
ac fields and ~ 10 /xV/cm Hz"^/^ and nT Hz~^/^ for dc 
fields. The sensitivity of the technique can be further en- 
hanced by employing the spin-echo pulse sequence |26j . 
which can be used e.g. to characterize the field of a mi- 
crowave stripline or map out the spatial distribution of 
electron spins. Finally, the method proposed here can be 
used for detecting weak dc and ac fields in the presence 
of high backround dc magnetic and electric field. 

The work was supported by NSERC of Canada and 
the NSF grant to ITAMP at Harvard University and the 
Smithsonian Astrophysical Observatory. 



states 1 2) and |3) become nearly degenerate at magnetic 
fields near B = B* , the magnitude of c?off is significatly 
enhanced near the avoided crossing depicted in Figure 1 
(b). The sensitivity to a dc electric field is given by 

100d,s{BQ,Eo)VnV^V^' 

and equals 1.58 x 10^'"' V/cm Hz^^^^ for SrF molecules 
with the density 10^^ cm~'^ in a magnetic field of Bq = 
545 mT and electric field of Eq = 2.5 kV/cm. 

The range of detectable frequencies of electromagnetic 
fields can be extended by using paramagnetic molecules 
of higher spin multiplicity. For example, molecules 
offer a series of tunable transitions that can be used to 
probe ac electric fields in the same way as the transitions 
in the molecules described above. This is illustrated in 
Fig. [4] that shows the energy level structure of NH(^S^) 
molecules, and the corresponding sensitivities. A large 
value of the rotational constant of NH and a series of 
tunable dipole-allowed transitions allow for the possibil- 
ity to cover continuously a broad range of detectable ac 
fields in the THz frequency region, which is particularly 
interesting for a variety of practical applications [25]. 

The sensitivity to the magnetic component of 
an ac field can be calculated as rylf [T Hz~^/^] — 
7?f [V/cm Hz-i/2] X 3.336 • 10"^ The ratio of the mini- 
mal detectable magnetic fields in experiments with cold 
atoms [Ej and molecules is given by B^-^/B^°^ — 
2cd]^ / (-yS/us) (assuming the same interaction time r, the 
same number of particles and the same detection effi- 
ciency in both cases), which for typical molecules with 
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